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Abstract 

A careful investigation of 1/rrib power suppressed effects is crucial for enhancing 
predictive power of the QCD factorization approach to charmless B decays. 
It is instructive to systematically investigate the l/m& effects from soft and 
hard gluon exchanges, in addition to annihilation topology in the charmless B 
decays. In this work we try to give a systematic discussion on impact of such 
soft and hard exchanges on the penguin-dominated B — > Kit decays within 
the framework of the light-cone QCD sum rules (LCSR). For the weak phase 
7(= ImV^) ranging from 40° to 80°, we find that because of the annihilation 
and soft and hard exchanges a numerical increase of (20 — 30)% is expected 
in the branching ratios up to 0(1/ ml) corrections; the resultant soft and hard 
corrections are less important than the annihilation contributions and amount 
only to a level of 10%. Possible sources of uncertainties are discussed in some 
details. 
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I. INTRODUCTION 



A considerable progress has been achieved in theoretical study of B decays since the naive 
factorization ansatz [1] was proposed for the non-leptonic decays of heavy-mesons. With the 
QCD background one makes all effort to approach the physics of B decays, by developing 
various theoretical frameworks such as large energy effective theory (LEET)[2], three-scale 
perturbative QCD factorization theorem [3], generalized factorization [4], QCD factorization 
[5] and soft-collinear effective theory (SCET) [6, 7, 8]. Especially, more attentions are paid 
to the QCD factorization and SCET. Both of them are formulated from the first principle 
of QCD and provide a rigorous theoretical basis to the factorization for a large class of non- 
leptonic B decays. In comparison, the QCD factorization approach turns out to be much 
simpler in structure and more intuitive in physical picture, while SCET possesses stronger 
predictive powers and provides a more complete theoretical framework, in which, for instance, 
an elegant proof of the factorization for B decays into two light mesons and into Dtt is given 
in Ref. [7] and some relations can be established among different decay modes [8]. 

The QCD factorization formula has been proved to be a substantial improvement of the 
naive factorization assumption. It allows us to compute radiative corrections to factorizable 
amplitudes to all orders in the heavy quark limit m& — > oo, meanwhile leaving the l/m& power- 
suppressed corrections to be estimated with help of concrete models. The reasoning behind 
this theory is that as one works in the heavy quark limit, the soft gluons with momentum 
of order Aqcd decouple and so the interaction kernel responsible for the transitions can be 
calculated with the perturbative QCD (PQCD) in the case of the B decays into two light 
mesons. To be specific, the hadronic matrix elements for, say, B — > nn can be expressed as, 



where Oi are the concerned local four-quark operators in the weak effective Hamiltonian, 
ji t 2 are the bilinear quark currents and the other two terms correspond to the perturbative 
corrections and non-perturbative contributions respectively. Nevertheless, an existing problem 
in the approach is that while the power corrections in a s can be calculated in a systematical 
way, the l/m& power-suppressed effects cannot. Thus an accurate theoretical prediction on 
nonleptonic B decays still is a challenge. 

In the QCD factorization, there are a variety of sources of power-suppressed effects. Among 
them the annihilation topology, soft and hard exchanges and final state interaction (FSI) 
effects etc. are the main ones. The recent work of Mantry, Pirjol and Stewart [8] indicates that 
there exists an additional source of power suppressed contributions. It is obvious that including 
all the power-suppressed effects is so far almost impossible in any practical computation. But 
it is plausible to make an order of magnitude estimate of the overall power-suppressed effects 




(1) 



by calculating the effects due to annihilation topology and soft and hard-gluon exchanges. To 
this end, a better understanding or at least a reliable estimation of the order of magnitude 
for both effects is crucial to enhance the predictive power of the QCD factorization. Based 
on PQCD [9], contributions of annihilation to B — > -nix, Kir have been estimated, showing a 
subleading behavior in 1/rrib. Soft effects can be understood as processes where a background 
field gluon is exchanged between, for instance, an emitted meson and the other meson which 
picks up the spectator quark in the case of emission topology. Therefore, they can be viewed as 
a higher-twist effect. Power suppressed hard effects originate from hard gluon absorption by a 
spectator quark, which is only relevant to the penguin contractions of effective operators and 
the chromomagnetic dipole operator. There have already been some earlier attempts [10] to 
understand such soft effects, and recently a systematic discussion is given by Khodjamiriani 
[11]. In Ref.fll] the author suggests using the QCD light-cone sum rule(LCSR), which is 
originally developed in Ref.[12], to evaluate the non-factorizable corrections to B decays into 
two light mesons. Using the generalized LCSR technique [11, 13, 14] non-factorizable soft 
corrections to B — > ixix have been discussed in details. It is found that despite its numerical 
smallness, the 0(l/m b ) soft effect is at the same level as the corresponding 0(a s ) corrections 
to the factorizable amplitudes. This implies that such soft effects are indispensable for an 
accurate evaluation of nonfactorizable contributions to charmless B decays. 

Recently, the CLEO-II and-III, Belle and BaBar Collaborations reported their data on the 
branching ratios of the B — > Kn decays [15]. The experimental averages are[16] 

(21.8 ± 1.4) x 1(T 6 , 
(18.2 ±0.8) x 10~ 6 , 
(12.5 ± 1.0) x 10~ 6 , 

(11.7 ± 1.4) x 10" 6 . (2) 

Contrast to B — > nil decays, these decay modes are penguin-dominated, for the tree contribu- 
tions are CKM suppressed. As well known, all new physics effects can only manifest themselves 
via loops, so that these processes deserve a detailed investigation in search for new physics. In 
this work we will make an evaluation of the power-suppressed effects in the B — > Kir decays 
by investigating all the possible soft-gluon effects, as well as the power-suppressed hard effects 
associated with the penguin topology and chromomagnetic dipole operator in the framework 
of the LCSR. 

Our paper is organized as follows. In the following section, we present a systematic LCSR 
analysis of the soft contributions to B — > Kn. The calculation is performed to 0(l/m b ). In 
particular, the soft effects stemming from emission topology are investigated at some length. 
Also, the special role of the penguin topology and chromomagnetic dipole operator is discussed. 



B(B~ -> K°tj-) = 
B(B° -> K-7i + ) = 
B(B~ -> K-n°) = 
B(B° -> K°n°) = 



They contribute, in addition to a soft effect, a f /m fe -suppressed hard effect which is absent 
in the QCD factorization. We estimate both of them with help of the LCSR results for 
B — > -kit [14]. In section III, we give the formula for calculating the branching ratios at 
subleading orders in a s and l/m&, including all the estimated power-suppressed corrections. 
In Sec. IV, we present our numerical discussions along with all the necessary input parameters, 
and a somewhat detailed analysis on the main sources of uncertainties, which influence our 
numerical results more or less. The last section is devoted to a brief summary. 

II. LCSR FOR SOFT-GLUON EFFECTS IN B -> Kn 

We begin with the weak effective Hamiltonian H e jf for the AB = 1 transition as [17] 

p 10 

n ^ts = -m E \ p (cmom p + c 2 (/i)o 2 (^ + Ea(/i)o l (/i) + c 83 ( / i)o 8 ,( / i)) + / i .c, (3) 

V ^ p=u,c i=3 

where \ = V pb V* q (q = d, s), Cj(/i) and C 8g (fj) denote the Wilson coefficients, 01,2, Oj(« = 
3 — 6) (Oi(i = 7 — 10)) and Sg are the tree, QCD (electroweak) penguin and chromomagnetic 
dipole operators, respectively. For a completeness, we list the relevant effective operators 
below 

0\ = (pb) v ^ A (sp)v- A , 

01 = {Pabf3)v-A{Sf3Pa)v-A, 

q 

4 = (s a b/3) V - A ^2(9l3Qa)v-A, 

q 

5 = (sb) v _ A J2(qq)v+A, 

q 

Oq = (s a bp) V - A ^2(qpqa)v+A, 
q 
3 

07 = (sfyv-A^ -e q (qq) v+A , 

q A 

3 

08 = (S a bp)v-Aj2 nZqiWc^V+A, 

q 1 
3 

9 = (sb) v ^ A J2 ~e q {w)v-A, 

q 1 

3 

010 = (s a bf3) V - A J2 o e q(Wa)v-A, 

q 1 

°S9 = ^m b sa^(l + 75)G^> ( 4 ) 

where a and (3 are the color indices, q runs over u, d, s, c and b, G^ v = G a is the gluon field 
strength. 



To estimate the soft-exchange effects from emission topology in the B — > Kix decays, 
one would decompose the relevant 4-quark effective operators into a color singlet part and a 
color octet one. Then the soft corrections to the matrix elements of the color octet operators 
can be estimated by studying the interactions between an emitted quark- ant iquark pair and a 
background field gluon. In the first place, we concentrate ourselves on the case of B° — > K~n + 
and calculate such effects in terms of the LCSR method. The operators, which may induce 
soft-exchange effects, contain the tree operator 2 , QCD penguin operators 3i5 , and EW 
penguin operators Oj t g. The color-octet operators which we encounter in this work, are of two 
types of structure: (V — A)(V — A) and (S + P)(S — P). The former case has been taken 
into account for B — > tttt in Ref. [11, 13]. Here we would like to provide a detailed derivation 
of the soft contributions arising from this type of operators in the B° — > K~n + case. To be 
definite, we are going to calculate the soft correction to the matrix element of the operator 

A a A a 

O = (u—b) v ^ A (s—u) V - A . (5) 

Given that the K~ meson is produced as an emitted hadron, the vacuum-pion correlation 
function, as the beginning point of the LCSR calculation, is written as 

F a (p, q,k) = -J d 4 xe-^-«>* J d 4 ye^- k >y(0\T{jif (y)O(0)ji(x)}\7r(q)), (6) 

where j^P = u^^^s and = mj,ib^d are the interpolating fields for K~ and B° mesons 
respectively. Then we decompose the correlation function (7) with respect to the independent 
momenta into four invariant pieces: 

F a (p, q, k) = {p- k) a F + q a F x + k a F 2 + e apXp q p p x k p P 3 . (7) 

Here an unphysical 4-momentum k ^ is introduced as an auxiliary external momentum in 
the weak operator vertex. Thus the total momentum of the final state becomes P = p — k — q, 
which is independent of the momentum p — q in the B channel. The advantage of introducing 
k is to help avoiding artificial ambiguities in the dispersion relation for the B-meson channel. 
Of course, the unphysical k has to vanish automatically in the physical matrix elements. This 
can be guaranteed, as will be seen, by choosing kinematical regions in such a way that we let 
A; 2 = 0. 

Saturating the correlator (6) with a complete set of intermediate states of the K quantum 
numbers and utilizing the definition of the K meson decay constant < | ir/o^s \K~(p— k) >= 
ifi<:{p — k) a , it follows that only the invariant function F is relevant to our concern. Explicitly, 
the resultant hadronic expression for F reads, 

Ph((P — k) ,(p-q) ,P)- > k _ {p _ k)2 +] s0 ds s _ {p „ k)2 , (8) 



where Sq is the threshold parameter and the spectral function p ( h stands for the higher state 
contribution in the K channel and U vK ((p — q) 2 , P 2 ) is a correlator, 

n^((p-g) 2 ,P 2 ) =ij d 4 xe i ^*(K-(p-k)\T{O(0)ji B \x)}\n-(q)). (9) 

On the other hand, applying the Operator Product Expansion (OPE) to Eq.(6) F can be 
calculated in large space-like regions of both p — k and p — q. In terms of the quark-hadron 
duality and then making the Borel transformation (p — k) 2 — > M 2 , we have 

K 2 

lW(p-g) 2 ,P 2 ) = — / dslmF QCD (s,(p-q) 2 ,P 2 )e^. (10) 

txJ k Jo 

The correlator U wK ((p — q) 2 , P 2 ) is applicable for only large space-like P 2 , therefore one needs 
to make an analytic continuation of Eq.(10) from a large space-like region P 2 <C to a large 
time-like region p 2 = m 2 B for the realistic B° — > K~n + decay. We then have 

U nK ((p-q) 2 ,M 2 B ) = i J d'xe^^K'ip - k)7r + (-q)\T{O(0)ji B \x)}\0) 

= — / dslmF QCD (s,(p-q) 2 ,M 2 B )e^. (11) 

71 f K Jo 

Inserting further the intermediate states of B° quantum numbers in the correlator of Eq (11), 
the hadronic matrix element < K~~n + \0\B > can be extracted in the light of the standard 
procedure for the QCD sum rule calculation. Not giving any technical details, we end up with 
the following LCSR expression for the matrix element in question, 

(K-(p)n + (-q)\0\B(p-q)) = — -/ dse^ ds'p Q c D (s,s',m 2 B ) 

JkJb^b Jo Jm 2 b 

m^j — s f 

x e^ 7 ^, (12) 

where fs is the B decay constant, Sq and M' 2 are the threshold and Borel parameters in the B- 
channel and the QCD double spectral density pQCD(s,s',m 2 B ) = l/7i 2 Im s Im s >FQ C D(s,s',m 2 B ). 
We will use the notation Ai to denote such matrix element from now on. 

Now we expand the correlation function (6) near the light-cone x 2 ~ y 2 ~ (x — y) 2 ~ in 
order to get Pqcd( s , s',m 2 B ). The kinematical regions we choose are summarized as follows: 

q 2 = 0,p 2 = m 2 K ,k 2 = 0, (13) 

and 

(p - k) 2 ~ (p - q) 2 ~ P 2 -> -oo, (14) 

where the LCSR calculation is efficient and self-consistent. Moreover, we set all the light quark 
masses to be zero and neglect all the terms proportional to order of 1/m 2 which emerge in the 
calculations. As has been mentioned, the soft effect is due to the interactions between emitted 



quark- ant iquark pairs and a background field gluon, thus the underlying quark propagator 
would receive a correction term [18]: 

0) = ^TZ^k } dv{(l - v)^G^{vx) + va^(vx)^} (15) 
^ ' o 

where n is the space-time dimension. 

Using Eq.(15) and parameterizing the nonperturbative QCD effects with the three particle 
distribution functions of the pion, a straightforward calculation leads to the following light- 
cone QCD result 

Fqcd = F tw3 + F tw ±, (16) 

with the twist-3 contribution 

m b f 37r f 1 , f n _ <P3k(oh) 



"tw3 = . 2 I dv / Da i7 - — 



An 2 Jo J (rn 2 — (p — q) 2 (l — «i))(— P 2 va 3 — (p — k) 2 (l — v a 3 )) 
x [(2-v)(q.k) + 2(l-v)q.(p-k)]q.(p-k), (17) 

and the twist-4 contribution 

t? m bf*(f\ fn ~< \ \ (2v-3)(p-k)-g 

F tw4 = — / dv / Da i (p ± (a i )- r ^ — r^iw ; ^~ 

n 2 \Jo J 2[m( — [p — (1 — ai)q) 2 \ (p — k — va 3 q) z 

f 1 dv [ daida 3 ^i(ai a 3 ) - (p ■ q — vq ■ k)(p — k) -q 

Jo J ' [ml — (p — (1 — cti)q) 2 } 2 (p — k — va 3 q) 2 

+ C dv /rfa 3 ^ 2 (tt 3 ) - {p ■ q - vq ■ k){p - k) -q 

Jo J [ml — (p — a 3 q) 2 ] 2 (p — k — va 3 q) 2 

f 1 . 2 f j t / \ 1 [{p~k)-qf 

- dvv cfo 3 ^ 2 (a 3 ) — — — 

Jo J p ■ q[m( - (p - a 3 q) 2 \ (p - k - va 3 qp 

In Eqs.(17) and (18), Da^ = daidai2da 3 8(l — ot\ — a 2 — a 3 ); f n is the decay constant 
of the pion, f 3w indicates a nonperturbative parameter defined by the matrix element < 
0\do' f j, u '~f5gsG a j3u\7i + >; ip 3n (ai) is a pionic twist-3 distribution function, whereas <f±(a{) is 
of twist-4 and it, together with another twist-4 distribution amplitude <p\\(ai), defines the 
functions ^i )2 as 

pu 

= J o dr] ((p±(m 1 -rj-v,v) + (p\\(m 1 -rf-v,v)), 

* 2 («) = / o " ^ jf ^ rfe (yM£, i - £ - v, v) + m(Z> 1 - £ - ^> ^) )■ ( 19 ) 



Readers are advised to refer to Ref.[19] for the definitions of various wavefunctions involved 
here and hereafter. 

A simple manipulation can make Eq.(16) expressed in a preferred form of dispersion in- 
tegral, from which extraction of p® CD (s, s' ,m 2 B ) is straightforward. Substituting the yielded 



ft CD (s, s',m B ) into Eq.(12), we derive the final LCSR results for A 1 as 

/ I fS K m%-s \ ( m 2 ,1 fa ml ml 



B \4n 2 f K Jo J \2f B m% J u * u 

m b f 3n f u dv f u dv. 



x 



u 



/ V37r(l - U,U ~ V,v) + fa / [3</?j_(l - U, U - V, v) 

JO V JO V 

(20) 



ml _ \ ^{l-^v) +f ( m 2 b n \ ^ 2 {u) 



\uM' 2 J u J Jn \uM' 2 ) u 2 

where Uq = m%/sQ. It is noted that for the twist-3 part we obtain the same result as that 

in B — > 7T7T case [11], whereas the obtained twist-4 parts are not quite the same (in contrast 

with the corresponding term (-^ - ^ - hi Ref. [11], our result is (^ - 2)*^). 

Numerically, however, the two forms result in close numbers. 

Applying the same procedure to the (S + P)(S — P) operators, one can notice that they 

do not result in any soft contributions to the amplitudes at all. 

Now let us turn to a discussion on the other three decay modes. The case of B~ — > K n~ 

is simple, where only the (V — A)(V — A) operators 03,9 are concerned so that the LCSR 

result (20) may apply directly. The situations for the B~ — > K~tt° and B° — > K°tt decays 

are a bit complicated, in which both K mesons may either be an emitted hadron or include 

a spectator quark (antiquark). It is obvious that the LCSR result (20) holds for the emission 

case where only the (V — A)(V — A) operators contribute to the decay amplitudes (note 

that there is an additional factor of l/\/2 (— l/y/2) for the operator (u^b)v^A(s^u)v~A 

((d^b)v-A(s^-d) V -A)- For the latter case, however, we have to modify the correlation function 

(6) with a necessary replacement, and besides the (V — A)(V — A) and (S+P)(S — P) operators 

we have to deal with the operators of (V — A)(V + A) structure. Omitting the concrete 

derivations to save space, we only present a simple summary of our results. The operator 

(s^-b)v~A(u^-u)v~A ((s^b)v-A(d^Yd)v-A) provides the relevant matrix elements associated 

with the soft contribution l/\f2A 2 (—1/V2A 2 ), 

... . *-o , „ s \ ( ml f 1 du nX-^k 
Ai = %m B I , - „ / dse m 7 | | — — I — £ M uM 



4tt 2 /. Jo J \2f B m% 



u 



X 



m b f 3K [ u dv r u dv 

/ — (p3k{l-u,u-v,v) + J K — [3if±{l - u,u- v,v) 
Jo v Jo V 

(21) 



u 

. m 2 b , ^i(l -u,v) ( m 2 c \ ty 2 (u 

•uM' 2 } u l+lK [uM' 2 J v? 



The (S + P)(S — P) operators make a vanishing contribution, as in the K emission case. It 
is interesting to notice that the (V — A)(V + A) operators have an equal matrix element to 
those of the corresponding (V — A)(V — A) operators. 

Expanding the LCSR results for A\ and A 2 in l/m& and then comparing them with the 
corresponding factorizable amplitudes, we observe that the soft effects of the emission topology 
are typically of order l/m b . 



Finally, the soft contributions of the emission topology to the B — > Kn decay amplitudes 
can be parametrized in terms of the resultant A 1 and A 2 as the following, 

M <P*\B- -> K\-) = -^V tb V*(2C 3 - C 9 )A U 

M^\B-^K-n°) = G F {V vb V^(C 1 A 2 + C 2 A 1 )-V fb V*[(C 3 + C Q )A 1 

-1(C S A 2 + C 10 )A 2 ]}, (22) 

MffB^W) = g f {k 6 K;Ci^2 + w;[(c 3 - ^C 9 )A! 

-l(C 8 + C 10 )A 2 ]}. (23) 

which evidently respect the isospin symmetry. 

Besides the soft contributions due to the emission topology, the two-body B decays, gen- 
erally speaking, receive the power-suppressed corrections from the chromomagnetic diploe 
operator 8g and penguin topology [14] . The relevant contributions contain a soft and a hard 
parts. The former is owing to a soft gluon which is emitted off either from the 8g vertex or 
from a quark loop of the penguin contraction and finally immerses into the meson which ab- 
sorbs a spectator quark. The latter is of two different origins. One is the hard gluon exchange 
between the 8g or penguin vertices and a spectator quark(antiquark). Another is related to 
the factorizable quark condensate contributions with QCD radiative correction being involved. 
The LCSR approach also is suitable for a quantitative study on these corrections. In fact, a 
detailed discussion has been made on their influences on B — > nn in the LCSR framework 
[14]. The behaviors of such contributions in the heavy quark limit m h — > oo comply with the 
following power counting: (1) The penguin diagrams make no contributions up to l/mf order. 
(2) The 8g operator supplies a hard effect of 0(l/m b ) and a soft effect of 0(l/ml). Despite 
being formally suppressed by OiAjmfy compared with the leading-order factorizable ampli- 
tudes, the soft correction from 8g is free of a s suppression and numerically comparable with 
the 0(a s ) hard part. (3) The part with quark condensate has a chiral enhancement factor r£ 
via the PCAC relation and it suffers formally from 0(a s ) and 0{l/m b ) double suppression 
but numerically turns out to be a large effect. However, it has a counterpart in the QCD 
factorization and thus is not included in our calculation to avoid double counting. 

It is reasonable to assume that the same power counting holds for B — > Kn. Furthermore, 
since no strange quark appears as a spectator in the concerned case, the emitted soft gluons 
can only combine with a quark pair to form a three-particle Fock state of the pion, we may 



directly use the corresponding LCSR results [14] by simply replacing the relevant parameters 
to achieve an estimate of the soft and hard effects due to the 8g operator and penguin 
topology in the B — > Kn case. As a consequence, the resultant corrections of Sg to the decay 
amplitudes can be expressed, at the subleading order in l/m^ as follows, 

= V2M ( ^\B- - K-n°) 

= -^V tb V: s C,Mf S9) + ^ a) )- (24) 

Here A { ° 8s) and A { ° 8s) express the hard and soft parts of the matrix element < K°n \0$g\B >, 
and have the following LCSR results: 




where both distribution amplitudes (fi p (u) and (p a (u) are of twist-3 and are used to describe 
the pionic valence Fock state, while y?j_(w) is a twist-4 three-particle wavefunction in analogy 
to &±{a,i) and (p\\(ai). 

III. DECAY AMPLITUDES AND BRANCHING RATIOS WITH 
SUBLEADING CORRECTIONS 

The LCSR results for the 0(l/m b ) soft and hard corrections to the B — > Kir decay 
amplitudes may serve as an order of magnitude estimate of the overall power-suppressed effects. 
We add them, together with the annihilation contributions M a , to the QCD factorization 
results Mf, to get a decay amplitude with the subleading power corrections in both a s and 
l/rrib. In Ref. [9], the B — > Kn decay amplitudes have been computed by including the 0(a s ) 
corrections, and the annihilation effects have also been estimated in PQCD. Their results, 
which will be used for our upcoming numerical discussion, are: 

M f (B~ - K\-) = £ V pb V; s [K - loSo) + rf K " \al)\A« K , 



Mf (B-^K-n°) = i=([K b K>i+ E V pb V; s (a p 4 + a p 10 ) 

V ^ p=u,c 



+ E VpbV p * s r«(a p 6 + a p 8 )}A wK 



p=u,c 



+ [V ub V: s a 2 + J2 VpbV p * s -(-a 7 + a 9 )]A K7T ), 



p=u,c 



M f (B° -> K-n+) = [KbK>i + E + o? ) 

P=M,C 

+ E K + At*, 



p=u,c 



and 



Mf(B° -> ifV) = Mf(B~ -> fsTvr ) - -^=M/(S~ -> A^vr") 

- ^,(5° (27) 



M a (B" - A^tf ) = ^M a {B~ - X°7T-), 

M a (B° - K-7T+) = (KfeKs + W)(6s - 

M a (5° -> AV) = -^M a (5° -> K-7T+). (28) 

v2 

In Eq.(26), and af are the QCD modified effective coefficients, and the low energy effects 
for the heavy to light transitions are included in the parameters A^^Ak^) with an obvious 
dependence on the B — > ir(B — > A) form factor F B ^ 7T (m 2 K )(F B ^ K (ml)), the parameter r£ 
is a so called chiral enhancement factor which is related to the running quark masses. As the 
parameters existing in Eq.(27) are concerned, b 2 and b 3 (b§ w ) are the parameters related to the 
tree operator O p and QCD (electroweak) penguin operators and embody the QCD dynamics 
in the annihilation processes, while B n x is closely related to the decay constants. The explicit 
definitions of all these quantities are given in Ref. [9] and we do not repeat them here. 

At present, we can write down the decay amplitudes of B — * Kir with the sub leading 
corrections in both a s and 1/m,},, 



M(B — > Kir) = Mf(B — > Kir) + M a (B — > Kir) 

AOs,. 

l h+s 



+ M^\B -> Ktv) + M { ° + 1 9 \b~ -> Att). (29) 



Using the resultant decay amplitudes, it is straightforward to calculate the branching ratios 
of the B — ■> Kn decays, which are given by 



B(B -> Kn) = ^\M(B -> 7Cr)| 2 ^, 



(30) 

87r m| 

where P is the cm momentum of the outgoing mesons in the center of mass frame of B meson, 

101 [( m s ~~ + m 7r) 2 )(m 2 B - (m K - m^) 2 )]^ 



2m i 



(31) 



and tb is the B lifetimes. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



Now we are in a position to make a numerical analysis and then to see how the subleading 
effects in l/m b impact the results of the QCD factorization. The set of important parame- 
ters of the present concern contains the QCD modified effective coefficients, quark masses, 
form factors, decay constants and distribution amplitudes. For the QCD modified effective 
coefficients and current quark masses, we adopt the numerical values presented in Ref. [9]. 
There have been a number of model-dependent estimates for the form factors F B ~* n ' K (q 2 ) 
and decay constants fs in the literature. Consistently we pick up, as inputs, the sum rule 
evaluations [20]: F B ^ n (0) = 0.28±0.06, F B ^ K (0) = 0.32±0.05 and f B = 180±30Me\/ with 
mj, = 4.7 ±0.1 GeV. The decay constants of 7r and K mesons are taken as f n = 132 MeV and 
fx = 160 MeV, respectively. As for the distribution amplitudes for the pion and K meson, we 
don't take into account SU(3) breaking effect for consistency, and adopt the following forms 
[19]: 

<P*,k(») = M 1 -«)[! + a 2 (fi)C 3 2 /2 (2u - 1)], 
<Pp(u,fj,) = 1 + 30^C 2 1/2 (2w - 1) - 3— ^C\'\ 



(Pcr(u,[x) = 6u(l - U) 

V^37r(«i) = SQ0aia 2 al 



1 + 5 



3tt 



^371-^3/2 
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)Ct'\2u - 1) 



1 + ^(7a 3 - 3) 



<P±(<Xi) = 30<5 2 (ai - a 2 )al[- + 2e 7r (l - 2a 3 )], 



<p\\ (ai) = -1205 2 aia 2 «3 
<P±(<Xi) = 30(^(1 - a 3 )aj 



- + e T (l - 3a 3 ) 



- + 26,(1 - 2a 3 ) 



(32) 



where the the leading twist- 2 distribution amplitudes (p n ,K(fJ>) enter in calculation of the QCD 
factorization and annihilation contributions. The various coefficients in above equation have 
been determined [19] at fi b = Jm 2 B - m\ ~ 2.4GeV: / 37r = 0.0026GeV 2 , u 3n = -2.18, 



fi w (= m 2 /(m u + m d )) = 2.02GeV, 5 2 = 0.17GeV 2 , e = 0.36. The scale-dependence of them is 
achievable by use of the renormalization group equations [14, 19]. The parameters inherent in 
the sum rules can be fixed via the two-point sum rules at: Sq = 0.7GeV 2 and M 2 = 0.5 — 1.2 
GeV 2 [11], sff = 1.2 GeV 2 and M 2 = 1.5 — 3 GeV 2 [21]; = 35 ± 2 GeV 2 and M' 2 = 8 — 12 
GeV 2 [20], corresponding to the pion, K and B channels respectively. The 5-lifetimes are 
measured experimentally as [22]: rgo = 1.542 x 10~ 12 ps and t b - = 1.674 x 10~ 12 ps. In 
addition, to explicitly investigate the dependence of the resultant branching ratios on the 
weak phase 7 = lmV* b , it is convenient to employ the following parametrization for the CKM 
matrix elements [22]: 

V ub V: s = (I - ^)\V cb \R b e-*\ 

v cb v; s = (1 - ^)\v cb \, 

V tb V* = -V ub V: s - V cb V; s . (33) 

where |V cb | = 0.0395 ± 0.0017, A = |V^| = 0.2196, and R b = ttVI^I with \V ub /V cb \ = 
0.085 ± 0.02. 

Before an explicit numerical analysis is made, we would like to simply judge the pattern for 
the rates of B — > Kir by analyzing the results of naive factorization. The following hierarchy 
among the rates should be expected: 

T(B~ -> K\-) > T(B° -> K-7i+) > T(B- -> R-ji ) > T(B° -> (34) 

The reasons are the following: (1) The penguin contributions dominate in these decays, since 
the tree ones are less important by the CKM suppression. (2) It is generally desirable that the 
branching ratios for B~ — > K°ir~ and B° — > K~ir + must be very close due to the smallness 
of the electroweak penguin effect, the slight difference between them arises from a destructive 
interference between the tree and QCD penguin contributions to B° — > K~n + . (3) The 7r° 
wavefunction makes the first two modes in Eq.(33) larger than the other two, and T(B~ — > 
K~ir°) nearly half as large as T(B — > Kn^). (4) Despite the fact that the electroweak 
penguin has a negligible effect on the first two modes, it plays a non-negligible role to the 
last two. With a moderate electroweak penguin contribution, the interference would become 
constructive between the electroweak and QCD penguins in B~ — > K~n°, in contrast to the 
B° — > K n° case where the interference is destructive between the two penguins. Qualitatively 
the pattern is in agreement with the experimental observations. 

Although the numerical calculation can be done using the inputs given above, we should 
note that the uncertainties resulting from such calculations are difficult to be quantitatively 
evaluated because of our poor knowledge about the inputs, especially the higher-twist distribu- 



tion amplitudes. However, the yielded results are adequate to serve as an order of magnitude 
estimate of the effects in question. 

we present our predictions based on the QCD factorization formula B^(B — > Kir) for 
each mode, respectively in Fig. 1-4, together with those with the annihilation effects included 
B^ +a \B — > Kir) and the results all the estimated subleading corrections in a s and l/m b 
BW(B -> Kir). The yielded B^(B -> Kir) and B^ +a \B -> Ktt) are plotted for the four 
different modes, respectively in Fig.5 and Fig. 6. It should be understood that only the central 
values, which are obtained for fi = fi b and 7 = 60° ± 20°, are shown there. Working within 
the QCD factorization framework, one observes from Fig.5: (1) The yielded branching ratios 
answer to the pattern 

B(f\K°7r-) > BV\K-ir + ) > B^(K-tt°) > B (f \K°n°). (35) 

(2) The branching ratios of both K°tt~ and B — > K°ir Q are less sensitive to the change of 7 
than the other two, as expected. 

As the l/rrib power-suppressed effects are involved from the annihilations and soft and hard 
gluons, the branching rations get an evident modification as shown in Fig.6. Letting \i vary 
between n b /2 and 2/i b , we find that the l/m b power-suppressed effects from the annihilations 
and soft and hard gluons can enhance the branching ratios by (20 — 30)%, depending on the 
decay modes and 7. 

Now, for taking a closer look at the individual roles of the soft and hard exchange and anni- 
hilation effects in the subleading contributions in l/m b , we estimate the two ratios B^ nV> /B^ +a ^ 
and £>(^ +a ) / B^ which can make sense about the involved physics. It is found that for 
/i ranging from /i{,/2 to 2fi b the effects of annihilation topology can make the branching 
ratios B^(B — > Kir) increase by (20 — 30)%, whereas the soft and hard effects modify 
#(/+»)(£ -> Kir) at a level of about 10%. 

As emphasized, most of the inputs suffer from theoretical uncertainties, which can affect 
the numerical results and should be carefully examined. Obviously, these uncertainties are 
related to the non-perturbative QCD, about which a solid knowledge is absent at present. 
Thus we are not so ambitious to make a complete quantitative estimate of them, instead, 
we just list the main sources of uncertainties which we can conjecture, and discuss how they 
influence the present results. By our observation, the most important sources of uncertainties 
are the following: (1) Distribution amplitudes. The models, which are here adopted for the 
distribution amplitudes of light mesons, are based on an expansion in conformal spin. For 
the leading twist-2 distribution amplitudes, the asymptotic forms are known precisely and the 
non-asymptotic corrections, in spite of their model-dependence, are also believed to be under 
control. In contrast, little is known about the uncertainties in higher-twist wavefunctions of 
light-mesons and distribution amplitudes of B mesons. Therefore, this could have a significant 



impact on the reliability of theoretical prediction. For example, the higher-twist wavefunctions 
can greatly affect the accuracy of the l/m b correction parts and the hard spectator scattering 
contributions, which are proportional to 1/Aqqd- When we carry out the calculation by using 
the extensively adopted form of the B meson wavefunctions, a considerable uncertainty may 
exist. (2)"Chirally enhanced" trems. A term proportional to the parameter appears in 
the calculations of the matrix elements of the (S — P)(S + P) operators. The part is called 
the "chiral enhancement" term, being formally l/m5— suppressed, but numerically large. Be- 
cause the factor is sensitive to the current quark masses whose precise values are not known, 
their variations may cause a considerable uncertainty in the numerical results. (3) Form fac- 
tors. One may believe that long-distance effects dominate the hadronic matrix elements of 
the heavy-to-light transitions from a naive power counting. If it is true, the LCSR results for 
the form factors, which are used in our calculations, should be relatively reliable. Neverthe- 
less, there are other viewpoints contrary to it [3] , that is, the short-distance contributions are 
predominant over long-distance ones so that PQCD is applicable in this case. A better under- 
standing of the uncertainties due to the form factors asks for a clarification of the transition 
mechanism. (4) Decay constants. In QCD sum rule calculations the decay constants for the 
B mesons are sensitive to the b quark mass, as we know. In contrast, the decay constants 
of the light pseudoscalar mesons have been experimentally measured to high accuracy, the 
uncertainties from this part is the least. (5) Final state interactions (F SI) . Besides those ef- 
fects which have already been estimated, the FSI effects may play a non-negligible role in the 
evaluation of decay amplitudes. These effects are usually regarded as being power suppressed 
in l/mb owing to the large energy released in the B decays. Unfortunately, so far, we lack 
a convincing quantitative evaluation on them. The B — > Kn modes can receive the final 
state rescattering contributions from a number of intermediate states. The FSI contributions 
from the B — > D S D followed by D S D — > Krr are anticipated to be especially important. The 
reason is obvious that the B — > D S D modes are CKM favored and thus can provide a large 
branching ratio; on the other hand, they are typically dominated by long-distance dynamics 
and violate the QCD factorization. (6) Other sources of power suppression. In SCET a new 
source of power suppressed contributions has been identified in study on the color suppressed 
process B° — > D n°. It is necessary to explore the corresponding impacts on B — > Kn. In 
addition, there might be other new sources of power suppressed contributions awaiting to 
be dug out and their effects should be estimated. (7) Theoretical approach. The numerical 
results presented here are based on a combination of the QCD factorization formula, PQCD 
approach and LCSR method. A possible inconsistency might manifest as the calculations are 
done by using all the three different approaches. Of course, to avoid this ambiguity is to 
do calculation within the same theoretical framework. The LCSR approach provides such a 



possibility. Practical manipulation, however, would be considerably difficult in view of the 
existing complications in technique. 

Once all these uncertain factors are taken into account and clarified, either theoretically 
or phenomenologically, with our new knowledge, the numerical results presented here can and 
should eventually be updated. 

V. SUMMARY 

We obtain the LCSR results for the l/m b suppressed soft and hard corrections to the 
B — > Kir decays. Then combining them with the previously obtained results of QCD fac- 
torization and estimates of annihilation topology, the branching ratios are calculated and the 
numerical impacts of such effects, as an order of magnitude estimate of the overall power- 
suppressed effects, are investigated to subleading order in l/m b . Also, the important sources 
of uncertainties are discussed in some details. 

Subject to a quantitative estimate of the existing uncertainties, for 7 = 60° ± 20° we have 
observed the following. (1) The l/mt, power suppressed effects resulting from the annihilations 
and soft and hard exchanges can make the branching ratios getting an increase of (20 — 30)% 
with respect to the results of the QCD factorization, depending on concrete decay modes. 
(2) The annihilation effect is predominate over the soft and hard ones, which modify the 
branching ratios by only 10% of the results which include the 0(a s ) radiative corrections and 
annihilation effects. 

Our present work helps to make sense about the contributions of the soft and hard-gluon 
exchange along with other power suppressed terms to the B — > Kn decays. No doubt, it still 
is too early to draw a decisive conclusion at present whether or not the theoretical estimates 
can accommodate the experimental data and the power-suppressed soft and hard effects are 
negligibly small in the B — > Kir decays. We have to await the improvement in experiment 
and progress in theoretical or phenomenological study on, amongst other things, the higher- 
twist wavefunctions as well as behaviors of FSI effects in the heavy quark limit. Certainly, 
a better understanding of the other sources of uncertainties is indispensable to arrive at a 
reliable conclusion. Anyhow a further investigation, whether theoretical or phenomenological, 
on the power-suppressed effects in charmless B decays may be needed. 
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Figure 1: B{B~ — > K°ir~) versus the weak phase 7. The solid, dotted and dashed lines 
correspond to B^(B~ -> K ^), B^ +a \B~ -> K°ir') and B {nl \B- -> K tt-), respectively. 
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Figure 3: B(B~ — > K~n°) versus the weak phase 7. The definitions of the lines are the same 
as in Figure 1. 
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Figure 4: B(B° — > K°ir°) versus the weak phase 7. The definitions of the lines are the same 
as in Figure 1. 
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Figure 5: B^(B — > Kir) versus the weak phase 7. The solid, dotted, dashed and dash-dotted 
lines correspond to B~ — > K°ir~, B° — > K~ir + , B~ — > i^~7r° and B° — > K 7r°, respectively. 
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